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Abstract: We study a strongly coupled Z2 lattice gauge theory with two flavors of quarks, invari- 
ant under an exact SU(2) x SU(2) x Ua(1) x 11^(1) symmetry which is the same as QCD with two 
flavors of quarks without an anomaly. The model also contains a coupling that can be used to break 
the Ua(1) symmetry and thus mimic the QCD anomaly. At low temperatures T and small baryon 
chemical potential fi b the model contains massless pions and massive bosonic baryons similar to 
QCD with an even number of colors. In this work we study the T — fiB phase diagram of the 
model and show that it contains three phases : (1) A chirally broken phase at low T and fis, (2) 
a chirally symmetric baryon superfluid phase at low T and high fis, and (3) a symmetric phase at 
high r. We find that the nature of the finite temperature chiral phase transition and in particular 
the location of the tricritical point that seperates the first order line from the second order line is 
affected significantly by the anomaly. 
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1. Introduction 



Understanding the QCD phase diagram in the temperature T and baryon chemical potential jiB 
plane is important from both phenomenological and experimental point of view [|l], ^ ^. Based 
on very general arguments it is possible to predict the important phases and the rough structure of 
the phase diagram [Q, |5|]. However, for a quantitative understanding of the phase diagram and to 
determine the location of critical points, careful non-perturbative calculations within the framework 
of QCD are necessary. While the finite temperature phase transition at = can be understood 
well ||6|, ^ ||, ^, 10], a reliable first principles calculation at ^ is impossible today due to 
the notorious sign problem. Despite this difficulty, many exploratory studies of the phase diagram 
have emerged over the past decade [11, 12, 13|, 0]. Since the systematic errors in all these studies 
are far from control, it is fair to say that important questions about the QCD phase diagram remain 



unanswered. For example, where is the critical point [15]? Can nuclear matter form a crystalline 



state []16|]? Can there be a state of matter where chiral symmetry is restored while the quarks 
remain confined (a quarkyonic phase) [|l7|]? Until a solution to the QCD sign problem is found 
these questions may not be answered satisfactorily. 

Given the difficulty of uncovering the QCD phase diagram from first principles, many studies 
have focused on model studies. These include random matrix models [18], 19, 20], NJL models [21, 
22, 23, 1^], PNJL and related models [25, 26, 27, 28], sigma models [29, 30] and strong coupling 
lattice models [0, ^^l- Some have also focused on two color lattice QCD [||, |4l, |5[ |3|, |7[ |38|] 



since it does not suffer from the sign problem and still has an interesting phase diagram [39]. Effects 



of the anomaly on the phase diagram have also been considered in model calculations [^, 41 , 42, 
43]. While model studies cannot give quantitative information about the QCD phase diagram, 
they are still useful to understand the nature of possible phases and critical points. Unfortunately, 
if one wishes to study these models from first principles using Monte Carlo methods, one again 



encounters sign problems. Hence most studies so far involve the mean field approximation. It 
would be useful to go beyond this approximation at least within models calculations. Recently, 
in a pioneering study, the T — phase diagram of the strong coupling limit of lattice QCD 
with staggered fermions was computed Q]. Some of the features of how nuclei emerge in lattice 
QCD was illustrated. Although the nuclear binding energy was found to be large compared to the 
real world, such ab initio calculations of nuclear matter from model field theories could teach us 
valuable lessons. 

In this work we continue this trend to perform ab initio calculations in model field theories 
and study a toy model of QCD with an even number of colors that contains massless pions and 
massive bosonic baryons. The T — phase diagram of QCD with an even number of colors is 
also expected to be complex and interesting. When the baryon chemical potential exceeds a critical 
value, baryons will condense to form a superfluid state. If this phase transition is first order, then 
the corresponding transition line can end on a second order critical point in the T — fiB plane just 
like in real QCD [15]. QCD with an even number of colors will also have exotic phases such 
as the quarkyonic phase and color superconductivity. When the number of colors is greater than 
two, baryons are naturally massive while the pions are massless. Two color QCD is special since 
baryons and pions are both massless due to extra symmetries that transform baryons into pions 
[p3|]. By breaking these symmetries baryons can be made massive while keeping pions massless 
in a lattice field theory. As far as we know, field theory models with massless pions and massive 
bosonic baryons have not been studied in the literature with Monte methods since they too suffer 
from sign problems in the conventional approach. 

World line formulations of lattice field theory offer new hope since they lead to new solu- 
tions to the sign problem [ 45 ] and allow us to accurately compute phase diagrams of lattice field 
theories with new algorithms p6| , |47| , ^ ^]. Recently, a model of two flavor QCD, invariant 
under SU{2) x SU{2) x Ua{^) x Ub{1), was constructed in the world line formulation [^, ^]. 
The model contains a coupling that can be used to break the anomalous Ua{1) symmetry and thus 
mimic the QCD anomaly. In this work we extend the model to include massive baryons. Since 
there is no sign problem we can study the T — ^b phase diagram accurately. We find that our 
model indeed contains a QCD-like tricritical point but its location strongly depends on the strength 
of the anomaly. We also observe that the model contains a confined but chirally symmetric baryon 
superfluid phase similar to the quarkyonic phase at high /i b and low T. 

Our paper is organized as follows. In section ^ we describe the model in detail along with its 
world line representation. We also discuss observables that we use to identify and distinguish the 
phases in the model. In section ^ we discuss our results of the phase diagram in the absence of 
the anomaly. In section ^ we discuss how the phase diagram changes in the presence of a large 
anomaly and how this change comes about as the anomaly is increased slowly. Section ^contains 
our conclusions. 



2. Model and Observables 



The model we study is a Z2 lattice gauge theory on a four dimensional hypercubic lattice whose 
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action is given by 



+-5 - ^2.1) 



where x denotes a lattice site on an x Lt hypercubic lattice and a = 1,2,3,4 denotes the 
direction. The fields ipx and ip^ are two component Grassmann valued fields given by 



-iAa; = ^ ^ J , '4^x= dxj (2.2) 

We choose ^q, = 1 for q = 1, 2, 3 and = T. The asymmetry factor T controls the temperature, 
the coupling c controls the strength of the anomaly, the parameter ms controls the mass of the 
baryon and the parameter /^^ is the baryon chemical potential. The parameter 5 is fixed based on 
convenience as discussed later. 



The most interesting feature of the action (2.1) is that it shares the global symmetries of two 



flavor QCD. When c = the action is invariant under a global SU{2) x SU{2) x Ub{1) x Ua{B) 
symmetry. Indeed it is easy to check that the action is invariant under the transformation 

V'x. ^ e^'^+^^^L^,,, V'x^ ^ e-^^^+^^'^iJVx., (2.3a) 
^ V^,^Lte-*^^-*^^ ^ i^xR^e^'^-^'^, (2.3b) 

where L,R € SU{2) and Xg and Xo refer to even and odd sites. When c / 0, the Ua{1) symmetry 
is explicitly broken and the action is invariant only under SU{2) x SU{2) x [/^(l). We interpret 
the C/a(1) symmetry as the anomalous axial symmetry of QCD and the parameter c will be used 
to change its strength of the anomaly [50, ^T|]. The C/b(1) is the baryon number symmetry and 



the chemical potential fiB couples to the baryon numbers as in QCD. In addition to the global 
symmetries the action is invariant under a Z2 gauge transformation ifj^ — )• —'4'x and ip^ —tpx at 
every site x. Thus, the baryons of the model are bosons made up of confined lid-diquarks. 

After integrating over the Grassmann variables, the partition function of the model can be writ- 
ten as a sum over world line configurations of neutral pions (uu, dd), charged pions (ud, du) and 
baryons {ud, du). Each configurations can contain five types of objects or fermion bags ^ ]. 
These are (1) instantons (isolated sites), (2) double-bonds (isolated bonds connecting neighboring 
sites), (3) neutral-pion loops (4) charged-pion loops, and (5) baryon loops. The loops are self- 
avoiding loops since they represent hard-core bosons. An example of a world line configuration 
with all the five types of fermion bags is shown in figure |l|. The weight of each configuration is 
the product of the weights of all the fermion bags in the configuration. Instantons have a weight c, 
double-bonds have a weight 2{ta)'^{l + e~^™^/2 — 6), and all the loops have a weight T"' where 
nt stands for the number of time-like bonds in the bag. The baryon loops are weighted with an 
additional factor exp(— m^n — ^bw) where n is the total number of bonds in the loop and w is the 
temporal winding of the loop. We will define 5 = (e"^™^ + 2 — ujb)) /2, so that the double-bonds 
have a weight uJoi^c 
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A variety of observables can be measured with our algorithm. However, for simplicity, we will 
mainly focus on the following three current-current susceptibilities: 



1. Chiral current susceptibility 



L3 

x,y 



2. Axial current susceptibility 



L3 



3. Baryon current susceptibility 

YB = ^Y.('^^^x)J^^y)). (2.6) 



L3 
x,y 



where J^i^) is the conserved chiral current, J^i^) is the conserved axial current and J^ix) 
is the conserved baryon current. Note that these currents are variables that are defined on the 
bonds connecting the site x with x + a. These currents can be measured easily in the world 
line configuration and Ref . [ 51 ] contains a discussion of how and are measured. can 
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similarly be measured by following the worldline of baryons. For example in figure ^ Ja{x) is 
zero on bonds that do not contain a baryon loop and on bonds which contain the baryon loop it 
takes value +1 or —1 depending on the direction of the baryon loop. The time component of the 
baryon current is the baryon number and its average over the lattice gives the baryon density ps- 

Using the current susceptibilities Yi defined above, it is easy to determine if the symmetry 
generated by the specific current ,i = C, A, B is spontaneously broken or not. One expects 

limy, = / Broken phase ^^.7) 

L-i-oo I Aexj){—aL) Symmetric phase 

In our analysis we will use this finite size scaling behavior to distinguish the phases. We will also 
need to distinguish between a first order and a second order phase transition between the various 
phases. This is accomplished by looking for the critical scaling relation 

LY.iL, T)^^fk m - T)Li/n'. (2.8) 

k 

which must hold close to a second order critical point Tc- In other words, if we plot LYi as a 
function of T for different values of L we expect all the lines to cross at a single point at a second 
order critical at T = Tc for large values of L. If this does not occur we claim the transition to be 
first order. In addition we also look for two state signals in the observables close Tc to confirm the 
first order nature of the transition. 



We update the world line configurations using the directed path algorithm discussed in 



Sip . We have tested the algorithm against exact results on a 2 x 2 lattice and the comparison is 
given in the appendix. Below we present results obtained on an x 4 lattice. We fix = 
0.1 and ujd = 10 throughout the study for convenience. For these parameters we find that the 
renormalized baryon mass Mb obtained through a baryon-baryon correlation function depends on 
T but is bounded in the range 0.8 > Mb > 0.5. 

3. Phase Diagram without Anomaly 

We begin with c = so that the action is invariant under the C/a(1) symmetry and the anomaly 
is absent. This is similar to QCD with an infinite number of colors. At small T and the 
SU{2) X SU{2) X Ua{1) symmetry is spontaneous broken into the diagonal SU (2) flavor symmetry 
leading to four Goldstone bosons. Since the baryons are massive the vacuum is free of them and 



the Ub{1) symmetry remains unbroken. Thus, as L increases, based on Eq. ( [2.7| ) we expect Yc 
and Ya will go to non-zero constants while will go to zero exponentially. At high temperatures, 
when all symmetries are restored, all three susceptibilities should go to zero exponentially. These 
expectations are consistent with our findings (see figure ^ where the /i^ = results are plotted). 



In order to identify the nature of the chiral phase transition we check if Eq. (2.8) holds. In 
figure |, we plot LYa,c{L) for L = 12,16,24,32 as function of T a.t fiB = 0, 0.2 and 0.4. 
Although there is a clear indication for a phase transition, it not consistent with second order 
scaling. In particular we do not see the lines at different values of L cross at a single point. If 
the transition was first order we also expect a two state signal at the transition point. In figure ^ we 
plot Ya and Yc at fiB = 0.0 as function of Monte Carlo time at T = 2.168 and L = 32. The two 
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fJ-B 


0.7 


0.8 


0.9 


1.0 




1.793(1) 


1.7879(6) 


1.7338(4) 


1.6485(3) 


xVdof 


4.39 


1.66 


1.39 


1.15 



Table 1: Results for Tc obtained from fitting to the XY scaling at c = 0. 



state signal is clearly visible confirming that the transition is indeed first order with Tc 2.168 at 
^iB = 0.0. 

For a fixed but small T, as fis increases we find that the baryon density pB jumps at a critical 
value of pb (which depends weakly on T). At the same critical value Yb also jumps to a non-zero 
value while Yc drops to zero. On the other hand Ya remains non-zero on both sides but also shows 
a small jump. As an illustration of these results, we show the data at T = 1.6 in figure ^. These 
results confirm that at large jj^b and small T the SU{2) x SU (2) symmetry is restored while [/^(l) 
and Ua{1) are both spontaneously broken. In other words the high ^b phase is a baryon superfluid. 
Since baryons in our model carry axial charge, it is not surprising that Ua{1) is also spontaneously 
broken along with [/^(l) in the superfluid phase. 

In figure ^ we plot the L dependence of Yi below the transition temperature (T = 1.6, left) and 
above the transition temperature (T = 1.9, right) at a fairly large chemical potential (pB = 0.8). 
Clearly, Yb and Ya are non zero in the superfluidity phase, but zero in the symmetric phase. Since 
the relevant symmetry now is Ua{1) x C/b(1), the phase transition can be second order in the 
universality class of the XY model [^]. Assuming the transition is second order we fit the Yb and 
Ya data close to the transition to Eq. (2^) at various values of > 0.7. Unlike the jj^b = case, 
now all the curves cross do cross at a point. We fix the critical exponent v = 0.6715 which is the 3d 
XY universal value [54]. The value of Tc at various fiB along with the /DOF is given in table |l]. 
The values of the other constants obtained in the fit for the Yb data at jj^b = 0.8 is /o = 0.308(5), 
/i = -0.149(2), /2 = 0.0161(1), /a = 0.0009(4) with fk = 0,k>4. The quality of the fit can 
be observed in figure where we plot LYb and LYa at ^b = 0.8 and fiB = 0.9. 

At fiB = 0.7 we find /DOF = 4.39 indicating the failure of the second order scaling. In 
figure P we plot LYb{L) as a function of T for different volumes at = 0.6. We observe two 
transitions, one atT ^ 1.58 from the chirally broken phase to the superfluid phase, and the other 




Figure 2: Plot of ¥,,, i = A,B,C asa function of L at c, = and T = 2.16 (left) and T = 2.19 (right). 
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1.92 1.94 1.96 1 



Figure 3: Plots of LYc (left) and LYa (right) as a function of T at different values of L. The data shown is 
for c = 0.0 and /is ~ 0.0 (top row) fiB ~ 0.2 (middle row) and jiB = 0.4 bottom row. 



at r = 1.75 from the superfluid phase to a symmetric phase. A second order finite size scaling 
analysis again gives us a large / DOF at both the transitions indicating that both are first order. 
The existence of a large fluctuations in the L = 32 data is due to tunneling between two metastable 
phases. 

Using the above information, in figure ^ on the left we sketch the entire phase diagram as a 
function of T and //^ in the absence of the anomaly. We find three phases, (1) A low T and low 

phase where SU (2) x SU (2) x Ua{^) is spontaneously broken into a flavor diagonal SU (2) 
symmetry. This is the chirally broken phase, (2) a low T large baryon superfluid superfluid 
phase driven by the spontaneous breaking of Ub{^) symmetry. Since baryons carry axial charge, 
Ua{^) symmetry also remains broken, and (3) a high temperature symmetric phase. For ^xb < 0.8, 
all the phase transitions are first order. There is a point M located roughly at T = 1.70(5) and 
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Monte Carlo Time 



Figure 4: Monte Carlo time evolution of Yq and Yji at L — 32, T = 2.168 and c,^b = 0. Tunneling 
between two meta-stable states which is characteristic of a first order phase transition is clearly seen. 




Figure 5: Plot of the Yc (top left), Yb (top right) Ya (bottom left) and the baryon density pB (bottom right) 
as a function of fiB st c — 0.0 and T = 1.6 for L — 32. 

Hb = 0.57(1) where all the three phases coexist. The first order finite temperature phase transition 
line from the superfluid phase to the symmetric phase ends at a tricritical point B after which the 
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Figure 6: Plots of Yc,Yb and Ya as a function of L at c = 0, /is = 0.8 and T = 1.6 (left) and T = 1.9 
(right). 



transition line turns second order in the XY universality class. 

Since the renormalized baryon mass changes with the paramter T we think it would be inter- 
esting to plot the phase diagram in terms of the dimensionless variable ^b/Mb{T). Further since 
T is also not the temperature but simply a parameter in the model, it would be better to use T/T^ 
where T^^ is the value of T at the phase transition when fis = 0- Hence, in figure ^ on the right we 
again plot the phase diagram in these scaled variables. Remarkably, we now find that superfluidity 
sets in when [ib/Mb ~ 1 which means that binding energy between nucleons, if any, is small in 
our model. 




Figure 7: Plot of LYb (left) and LY^ (right) at fiB — 0.8 (top row) and fis = 0.9 (bottom row) as a 



function of T. The solid lines are joint fits to Eq. (2.8) with ly = 0.6715 fixed to the 3d XY universal value. 
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Figure 8: Plot of LYg as a function of T at /i^ = 0.6. Two first transitions are seen: one from the chirally 
broken phase to the baryon superfluid phase at T 1.58 and another from from superfluid phase to the high 
temperature symmetric phase at T ^ 1.75. 

I 1-10 I 




Figure 9: Sketch of the phase diagram at c = as a function of T and (left) and as a function of rescaled 
variables T /Tc{^b — 0) and ^ib/Mb (right). Solid black lines first order lines and the dashed lines are 
second order lines. The point M is a three phase coexistence point and the point B is a tricritical point. 



4. Phase Diagram with Anomaly 



Next we set c = 0.3 so that the Ua{1) symmetry is explicitly broken by a large amount. Now, the 
action is only invariant under an SU{2) x SU (2) x Ub{1) symmetry. At small T and fiB we expect 
the SU (2) X SU{2) symmetry to be spontaneously broken to the diagonal SU{2) flavor group. As 
the temperature increases a phase transition from the chirally broken phase to the symmetric phase 
must occur at a critical temperature. If this phase transition is second order it would belong to the 
3d 0(4) universality class [50, 54 1. We have verified that our data fits well with this expectation 
for 0.0 < < 0.5 by fitting the Yc data to the scaling form given by Eq. (2.8). We fix = 0.745, 
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the 0(4) critical exponent, in the fits. The fitting results are tabulated in table ^. In figure [I^ we 
plot the data along with the fits at fiB = 0.2 and fiB = 0.5 in order to show the quality of the fits. 



fJ-B 
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X^/DOF 
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2.4568(4) 
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2.4078(6) 
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2.3200(5) 
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2.1958(4) 


1.20 


0.5 


2.0351(5) 
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Tc 


X^/DOF 


0.6 


2.0557(6) 


1.68 


0.7 


2.0811(6) 


1.70 
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2.0490(5) 


1.30 


0.9 


1.9715(5) 


1.20 


1.0 


1.8609(4) 


1.09 



Table 2: Fitting results of Yc for 0(4) scaling 
at c = 0.3. 



Table 3: Fitting results of Yg for XY scaling at 
c = 0.3. 



Similar to the c = 0.0 case, we again find that at low T as fis increases there is a first 
order phase transition from the chirally broken phase to a chirally symmetric baryon superfluid 
phase where Ub{1) is spontaneously broken. Again, the symmetry is restored at high temperatures 
through a second order phase transition which belongs to the 3d XY universality class for suffi- 
ciently high hb- If we repeat the scaling analysis of Eq. ( p^ ) on our Yb data with u = 0.6715 (the 
3d XY critical exponent) we again find excellent fits for /i^ > 0.6. The results are given in table ^. 
As an illustration of the quality of the fits we show the joint fit results at hb = 0.7 and hb = 0.9 
in figure [Tl] 

Tables |^ and || show the location of Tc as a function of obtained from the O (4) scaling and 
the XY scaling respectively. What happens in the region 0.5 < hb < 0.6 is unclear. Unlike the 
results in the previous section we are unable to find a region where the two second order transitions 
clearly become first order. Can the two different second order transitions meet at a multicritical 
point? This seems unlikely, but cannot be ruled out with our limited statistics. Most likely the 
transition lines becomes weakly first order before meeting. Based on this guess, in figure 12 we 



sketch the phase diagram of the model at c = 0.3. As compared to figure ^ the phase diagram in 
figure contains the extra tricritical point C. Through a more extensive calculation, the precise 
nature and location of the points M, B and C can be understood. 
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Figure 11: Plots of LYb as a function of T at fiB = 0.7 (left) and — 0.9 (right). The solid lines are joint 



fits to the scaling form given in Eq. (2.8) with ly — 0.6715 fixed to the 3d XY universal value. 



T 



2.6 
2.4 
2.2 
2.0 
1.8 
1.6 
1.4 



1.2 1 

0.0 



Symmetric 




Chirally 
Symmetric 
Baryon 
Superfluid 



0.8 



1.0 



1.10 
1.05 

1 .00 _ , ^(4) 
0.95 
0.90 



rO 0.85 
c 

0.80 
0.75 
0.70 
0.65 



Symmetric 



Brol<en 
Chiral Symmetry 



M 



C 



B XY scaling 

Chirally 
Symmetric 

Baryon 
Superfluid 



0.0 0.2 0.4 0.6 0.8 1.0 
^^B/MB 



1.2 1.4 1.6 



Figure 12: Sketch of the phase diagram at c = 0.3 as a function of T and iib (left) and as a function of 
rescaled variables T/Tc{^b — 0) and hb/Mb (right). Solid black lines represent first order Unes while 
dashed line represent second order lines. The point M is a three phase coexistence point and the points C,B 
are tricritical points. The precise nature and location of M, C and _B is a guess represented by a question 
mark. 



Having established the phase diagrams for zero anomaly and a large anomaly, we next focus 
on how the phase diagram changes from figure ^ in the absence of the anomaly to figure 12 in the 
presence of a large anomaly. Note that two phase diagrams are very similar except for the nature of 
the finite temperature transition from the chirally broken phase to the symmetric phase. Without the 
anomaly the transition is first order, while in the presence of a large anomaly the transition becomes 
second order in the 0(4) universality class. So how does the transition change from first order to 
second order as the anomaly is slowly increased? Different scenarios have been proposed based 
on an effective theory approach [pO|]. The conventional scenario is that the first order transition 
becomes weaker at /i^ = until it becomes second order above a critical value of the anomaly 
strength. Then, as the strength of the anomaly further increases the second order tricritical point 
C moves into the T — plane. However, an alternative exotic scenario suggests that the first 
order transition becomes weak somewhere in the middle of the first order curve. As the anomaly 
increases further the second order line appears in the middle of the first order line creating two 
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tricritical points (see figure 4, in [30]). Wiiat liappens in our model? Motivated by this question, 
we have studied the model at c = 0.005, 0.01, 0.02, 0.03 and 0.05. We have evidence that in our 
model the change in the phase diagram follows the conventional scenario. To demonstrate this we 
discuss the results at /x^ = 0.02 below. 

Figure |l3| shows that close to the finite temperature phase transition, Yc scales according to 
0(4) universality at /ub = 0.0 and 0.1. But this scaling breaks down at /is = 0.3 and 0.5 as can 
be seen in figure |l4[ In particular the large fluctuations in the L = 32 data is consistent with the 
values fluctuating between the two metastable phases. Figure |l^ shows the time histories of the Yc 
data at L = 32 for /i^ = 0.3 and fig = 0-5. These time histories clearly show the two state signal. 
Thus, we predict a tricritical point at roughly jj^b 0.25(5). On the high chemical potential side, 
the anomaly has little effect on the phase diagram once it is plotted in the scaled variables (compare 
the right figures ^and 13). It does change the location of the point B but only slightly. Thus, the 
phase diagram at c = 0.02 at large fiB should be very similar to the one at c = 0.0. Our data 
confirms this. 



5. Conclusions 

In this work we have computed the T — fiB phase diagram of a lattice model which is invariant 
under an SU{2) x SU{2) x C/a(1) x Ub{1) symmetry and contains massive bosonic baryons. 
We also studied the effects of breaking the Ua{1) symmetry. Due to the global symmetry and 
the nature of the baryons, our model may be considered as an interesting toy model of QCD with 
an even number of colors and may share some of the important phases and phase transitions with 
it. The term that breaks the Ua{1) symmetry is mapped to the anomaly in QCD. Based on the 
results presented in this work, we conclude that the phase diagram contains at least three important 
phases : (1) a low temperature chirally broken phase with massive baryons (2) A chirally symmetric 
baryon superfluid phase at a moderately high chemical potential and low temperatures and (3) a 
symmetric high temperature phase. A qualitative sketch of the phase diagram at a generic value of 
the anomaly strength is shown in figure [l^. We find that the location of the tricritical point C in the 
figure is strongly dependent on the strength of the anomaly. There is also a tricritical point on the 
baryon superfluid side (point B) and a point M where all the three phases coexist. If the points B 



0.6 




o0.3 




Figure 13: Plots of LYc as a function of T at /i^ =0 (left) and = 0.1 (right). The solid lines are joint 
fits to the 0(4) scaling form given in Eq. (|2.8|) with fj,B fixed to the value obtained from 3d 0(4) universality. 
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Figure 14: Plots of LYc as a function of T at — 0.3 (left) and 1.1 b = 0.5 (right). The data does not fit 
well to the 0(4) scaling form and is consistent with a first order transition. 




200 400 600 800 1000 
Monte Carlo Time 



Figure 15: Monte Carlo time history of yc data at /us = 0.3(left) and /xb = 0.5 (right). The two state signal 
is clearly visible. 



and M exist in real QCD the physics close to them could also be very interesting both theoretically 
and phenomenologically. 

In the presence of a quark mass the finite temperature 0(4) transition will turn into a cross 
over, while the tricritical point C will turn into a critical end point. The location of a similar critical 
end point in QCD is the main focus of many calculations today including one of the goals of the 
experimental program at RHIC [Q]. Among the many dynamical effects that affect its location, 
here we find that the strength of the anomaly is an important one. Given that the anomaly in lattice 
QCD calculations is strongly dependent on the lattice spacing, we believe it will be rather difficult 
to find the precise location of the critical point in real QCD without controlling all the systematic 
errors. On the other hand, the two other critical points B and M seem rather robust and may be 
easier to locate. 

The present work can be extended further in many directions. We believe that the simplicity of 
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Figure 16: Sketch of the phase diagram at a generic value of the anomaly strength. The first order critical 
lines are shown with solid critical lines are shown with a dashed line. The location of the critical point C 
in the diagram is strongly dependent on the strength of the anomaly, while the location of the critical lines 
(both first and second order) and other critical points are insensitive to it. 



the phase diagram in our model comes from the fact that baryons in our model were singlets under 
chiral transformations. One can build and study lattice field theory models where baryons transform 
non-trivially under the chiral symmetry. It would be interesting to see if these models contain 
a phase where both chiral symmetry and the baryon number symmetry are both spontaneously 
broken. Such a phase would make the phase diagram even more interesting. One can also compute 
the properties of nuclei and their physics within our toy model. It would also be interesting to 
understand the connection between the strong first order transition as a function of the chemical 
potential and the binding of the baryons. Finally we wish to emphasize that new methods of 
analysis to locate critical points in QCD, can first be tested in toy models such as the one studied 
here before being applied to QCD. 
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A. Comparison between exact results and the algorithm 

In order to confirm that the algorithm reproduces the results in the model accurately, we have 
obtained analytic expressions for the partition function, the axial current susceptibility (Ya), the 
baryon current susceptibility (Yb), and the chiral current susceptibility (Yc) an a 2 x 2 lattice. 
These are given below: 

Z = 4iojDf{l + T^) + 8a;D(2{l + T^} + ^-^^^ {I + cosh(2^B)}) + 16(1 + AT"" + T^) 
+16e-2™s{i + cosh(2^B)} + 4e-^™s{l + 4[1 + cosh(2/iB)] + cos\v^{2^lB)} 
+4.c^ud{1 + T^) + 4c^(2 + 2T2 + e'^^^jl + cosh(2^B)}) + (A.l) 



ZxYa = 8{ujDf + 8ujd{2 + + 8T^{2 + e'^^^) + 8T^[2 + e"^™^ cosh(2^B)] 

+4:C^ujD (A.2) 

ZxYb = 4e-^™^ {2ujD + 4 + c^) + Se"^""^ [l + T^ + cosh(2/iij)] (A.3) 

Z xYc = 16 + 8e-2™^ +8ujd + IQT^ + 4c^ (A.4) 

In tables ^ and ^, we compare the analytic results with the results obtained using our algorithm for 
different parameters. The agreement gives us confidence that our algorithm must be correct. 
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ruB 




Exact 


Algorithm 


Exact 


Algorithm 


Exact 


Algorithm 






Ya 


Yb 


Yc 








T = 


1.0, c = 0.0, UJD = 0.0 






0.0 


0.0 


0.28571429 


0.28573(2) 


0.23809524 


0.23808(2) 


0.23809524 


0.23808(2) 


0.0 


0.5 


0.27419502 


0.27420(2) 


0.23228892 


0.23228(2) 


0.20953053 


0.20951(2) 


0.0 


1.0 


0.22684802 


0.22684(2) 


0.20095862 


0.20096(2) 


0.12944702 


0.12947(1) 


0.2 


0.0 


0.28938325 


0.28939(2) 


0.15882791 


0.15883(1) 


0.27589442 


0.27589(2) 


0.2 


0.5 


0.27867066 


0.27869(2) 


0.15891388 


0.15892(1) 


0.25307451 


0.25305(2) 


0.2 


1.0 


0.23552474 


0.23552(2) 


0.15074712 


0.15076(1) 


0.17915523 


0.17915(2) 








T = 


1.0, c = 0.0, UJD = 1.0 






0.0 


0.0 


0.35714286 


0.35715(2) 


0.21428571 


0.21427(1) 


0.21428571 


0.21428(1) 


0.0 


0.5 


0.33570326 


0.33569(2) 


0.20833887 


0.20833(1) 


0.19104659 


0.19105(1) 


0.0 


1.0 


0.26374692 


0.26376(2) 


0.18108186 


0.18109(1) 


0.12399759 


0.12400(1) 


0.2 


0.0 


0.36826413 


0.36824(2) 


0.14435557 


0.14434(1) 


0.24369015 


0.24366(1) 


0.2 


0.5 


0.35029205 


0.35030(2) 


0.14320792 


0.14323(1) 


0.22537934 


0.22539(1) 


0.2 


1.0 


0.28633667 


0.28636(2) 


0.13426676 


0.13426(1) 


0.16550478 


0.16549(1) 








T = 


1.0, c = 0.2, UJD = 0.0 






0.0 


0.0 


0.28408822 


0.28413(2) 


0.23768714 


0.23769(2) 


0.23768714 


0.23767(2) 


0.0 


0.5 


0.27269729 


0.27272(2) 


0.23185364 


0.23186(2) 


0.20921955 


0.20922(2) 


0.0 


1.0 


0.22582231 


0.22582(1) 


0.20056541 


0.20057(2) 


0.12937716 


0.12939(1) 


0.2 


0.0 


0.28756535 


0.28756(2) 


0.15861716 


0.15861(1) 


0.27533532 


0.27534(2) 


0.2 


0.5 


0.27695539 


0.27695(2) 


0.15865772 


0.15866(1) 


0.25259387 


0.25258(2) 


0.2 


1.0 


0.23423048 


0.23423(2) 


0.15043018 


0.15042(1) 


0.17893372 


0.17892(2) 








T = 


1.0, c = 0.2, UJD = 1.0 






0.0 


0.0 


0.35582134 


0.35581(2) 


0.21377689 


0.21376(1) 


0.21377689 


0.21376(1) 


0.0 


0.5 


0.33451804 


0.33453(2) 


0.20784361 


0.20784(1) 


0.19064501 


0.19065(1) 


0.0 


1.0 


0.26298931 


0.26298(2) 


0.18069067 


0.18070(1) 


0.12385944 


0.12386(1) 


0.2 


0.0 


0.36680615 


0.36682(2) 


0.14402179 


0.14402(1) 


0.24301429 


0.24302(1) 


0.2 


0.5 


0.34894700 


0.34894(2) 


0.14286463 


0.14286(1) 


0.22479567 


0.22482(1) 


0.2 


1.0 


0.28538737 


0.28536(2) 


0.13393857 


0.13393(1) 


0.16520110 


0.16519(1) 



Table 4: Comparison of exact results with results from the algorithm at T = 1.0. 
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Algorithm 
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Algorithm 


Exact 


Algorithm 






Ya 


Yb 


Yc 








T = 


1.5, c = 0.0, UJD = 0.0 






0.0 


0.0 


0.24870466 


0.24870(2) 


0.13816926 


0.13816(1) 


0.13816926 


0.13816(1) 


0.0 


0.5 


0.22401091 


0.22400(2] 


0.13271409 


0.13272(1) 


0.11412102 


0.11413(1) 


0.0 


1.0 


0.15405340 


0.15406(1] 


0.10716914 


0.10715(1) 


0.05860532 


0.05859(1) 


0.2 


0.0 


0.25245829 


0.25246(2] 


0.08731264 


0.08730(1) 


0.16423607 


0.16422(1) 


0.2 


0.5 


0.23135853 


0.23135(2] 


0.08719658 


0.08719(1) 


0.14336794 


0.14337(1) 


0.2 


1.0 


0.16653699 


0.16655(2] 


0.07961861 


0.07962(1) 


0.08643966 


0.08644(1) 








T = 


1.5, c = 0.0, UJD = 1.0 






0.0 


0.0 


0.23178808 


0.23177(2] 


0.11258278 


0.11258(1) 


0.11258278 


0.11258(1) 


0.0 


0.5 


0.20874557 


0.20875(1] 


0.10839747 


0.10839(1) 


0.09477320 


0.09477(1) 


0.0 


1.0 


0.14533254 


0.14534(1] 


0.09017774 


0.09019(1) 


0.05209065 


0.05209(1) 


0.2 


0.0 


0.23366106 


0.23366(2] 


0.07128440 


0.07128(1) 


0.12993783 


0.12994(1) 


0.2 


0.5 


0.21443609 


0.21447(2] 


0.07078745 


0.07081(1) 


0.11495121 


0.11498(1) 


0.2 


1.0 


0.15675114 


0.15675(2] 


0.06521854 


0.06523(1) 


0.07324665 


0.07325(1) 








T = 


1.5, c = 0.2, UJD = 0.0 






0.0 


0.0 


0.24781350 


0.24784(2] 


0.13804130 


0.13804(1) 


0.13804130 


0.13804(1) 


0.0 


0.5 


0.22326474 


0.22326(2] 


0.13257534 


0.13258(1) 


0.11404420 


0.11404(1) 


0.0 


1.0 


0.15366975 


0.15366(1] 


0.10705814 


0.10704(1) 


0.05861527 


0.05862(1) 


0.2 


0.0 


0.25145743 


0.25147(2] 


0.08727235 


0.08727(1) 


0.16404125 


0.16404(1) 


0.2 


0.5 


0.23048223 


0.23050(2] 


0.08713335 


0.08713(1) 


0.14322330 


0.14322(1) 


0.2 


1.0 


0.16602244 


0.16601(1] 


0.07953373 


0.07952(1) 


0.08641295 


0.08640(1) 








T = 


1.5, c = 0.2, UD = 1.0 






0.0 


0.0 


0.23125599 


0.23124(1] 


0.11246010 


0.11245(1) 


0.11246010 


0.11246(1) 


0.0 


0.5 


0.20830747 


0.20832(1] 


0.10827684 


0.10828(1) 


0.09469566 


0.09470(1) 


0.0 


1.0 


0.14511317 


0.14513(1] 


0.09008803 


0.09009(1) 


0.05209047 


0.05209(1) 


0.2 


0.0 


0.23309400 


0.23310(1] 


0.07122714 


0.07123(1) 


0.12976315 


0.12976(1) 


0.2 


0.5 


0.21394944 


0.21395(1] 


0.07072219 


0.07072(1) 


0.11482041 


0.11482(1) 


0.2 


1.0 


0.15647860 


0.15648(1] 


0.06515061 


0.06515(1) 


0.07321454 


0.07321(1) 



Table 5: Comparison of exact results with results from the algorithm at T = 1.5. 
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